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Abstract 

A new mean field statistical mechanics model of two interacting groups of spins is 
introduced and the phase transition studied in terms of their relative size. A jump 
of the average magnetization is found for large values of the mutual interaction 
when the relative percentage of the two populations crosses a critical threshold. 
It is shown how the critical percentage depends on internal interactions and on 
the initial magnetizations. The model is interpreted as a prototype of resident- 
immigrant cultural interaction and conclusions from the social sciences perspectives 
are drawn. 

1 Introduction 

In the last few decades the statistical mechanics approach has seen an impressive expan- 
sion in fields as diverse as combinatorial optimization, finance, biology and others. Its 
success relies upon the fact that every problem of many interacting parts may be studied 
by its methods and techniques. Our plan in this work is to introduce a statistical mechan- 
ics model with the purpose of describing two homogeneous groups of individuals whose 
interaction is imitative within the same group but may be imitative or counter-imitative 
between different groups. Such a model may represent a first attempt toward the construc- 
tion of a statistical mechanics theory of the phenomenon of cultural interchange. When 
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two populations come in contact, like in the case of immigration but also in a more general 
context through the media, it is often seen how cultural traits are sometimes mixed, while 
some other times one of the two dominates and prevails on the other. Examples are found 
in linguistics, as well as in opinion forming dynamics (PP,[2],[3])- One interesting feature 
of those changes is that sometimes they appear to be smoothly varying with the relative 
proportions of the two groups, some other times the crossing of a critical value triggers a 
jump of the observed quantity [4J. What we are building here is the mean field theory of 
the two population problem i.e. we consider every individual interacting with every other 
with the same strength. In future works we plan to introduce a more realistic model by 
allowing randomness of the interaction strength like in the Sherrington-Kirkpatrick spin 
glass model ([S] [S]) and also a more structured network connection like, for instance, 
the one predicted by the small- world theory [7] [8]. We want to stress that, although 
our model is inspired by the Curie- Weiss theory of ferromagnetism, the problem we deal 
with here is quite different because we do not study the phase transition in terms of the 
temperature but in terms of the relative size of the two populations. 

The dictionary we plan to follow is easily explained by saying that a cultural trait is 
considered for simplicity as a dichotomic variable cij = ±1. The interaction between par- 
ticles is built up as a sum of pairwise interactions and plays the role of cultural interaction 
between two individuals i and j as described by a potential, or a cost function, which sim- 
ply reflects the will to "agree" or "disagree" among the two. The two attitudes of imitation 
or counter- imitation lie on a well established socio-psychological background [Sllinilll]; on 
the other hand they have also a robust mathematical-physical ground since they have been 
used to study many particles theory of ferromagnetic and anti-ferromagnetic interactions. 

The problem we have addressed with the help of an equilibrium statistical mechanics 
model is to establish whether -in the case of two populations placed in contact- there 
may be a phase transition in the average cultural trait from one of the two original 
cultural traits to the other. If so, for which value of the relative percentage of the two 
populations it happens. Moreover, we want to establish how the critical size depends 
on the original parameters in order to predict or potentially avoid unwanted dramatic 
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phenomena sometimes occurring in society. 

The parameters describing our system are ml and i.e. the magnetizations of the 
two populations prior to their interaction which represent the two culture legacies, the 
couplings Ji^i, J2,2 which measure the strength of the imitation within each group and Ji^2 
which measure the strength of the imitation or counter-imitation among the two groups. 
The phase transition is tuned by the parameter a = Ni/N which counts the percentage 
of immigrants, 1 — a = N2/N being the fraction of residents. 

Our results, explained in detail in section 3, show that when the mutual interaction 
between the two groups Ji 2 is small enough the transition from the resident to the immi- 
grant culture is smooth. But for large values of the interaction there is a critical value of 
the immigrant percentage Oc crossing which the system undergoes a sudden change from 
the resident to the foreign culture. 

We find moreover that high values of the culture legacy favour both the emergence of 
the immigrant culture (oc decreases with ml) and the persistence of the local culture 
increases with mg) as intuition would suggest. On the contrary, a high internal imitation 
(high coesion and low diversification) makes each culture weaker toward the other (etc 
increases with Ji 1 and decreases with J2,2)- This last result is rather counter-intuitive but 
not surprising since in social sciences it is often seen how a diversified culture dominates 
an opinionated one. From a technical point of view the dependence of ac on Ji i and J2,2 
is explained by the detailed balance between energy and entropy, the second being the 
leading term in those situations in which a phase transition occurs. 



see [12] for some results on the model. 

The symbol (jj represents the opinion of the i^^ individual of the total population. 



2 The Model 



The model we introduce is defined by the Hamiltonian 
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which can either take value ai = +1 or ai = —1. We consider only the case of symmetric 
matrices Jjj. The general case can be easily reduced to the symmetric one by standard 
methods. 

We divide the total population I into a partition of 2 subsets hD I2 — /, of A^i = 
and = I/2I with A^i + A^2 = N. Given two individuals and aj , their mutual 
interaction parameter J^j depends on their respective subsets, as specified by the matrix 

Ni N2 



( 



'12 



I v 



'12 



'22 



The interacting system is therefore described by three degrees of freedom: Jn and 
J22 tune the interactions within each of the two subsets, and J12 controls the interaction 
between two individuals belonging to different sets. We assume Jn > and J22 > 0, 
whereas J12 can be either positive or negative. 

Analogously, the field hi takes two values h\ and hi-, depending on the subset containing 
(Tj, as described by the following vector: 



N2 < 



ho 



( V / 

By introducing the magnetization of a subset / as 



mi 



and indicating by mi and m2 the magnetizations within the subsets Ii and I2 and by 
a — — the fraction of the first group on the whole, we may easily express the Hamiltonian 
per particle as 
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^ ^ = — - [JiiQ;^m^ + 2Ji2Q;(l — a)mim2 + ^22(1 — Q;)^^2] ~ hiarrii — /i2(l — a)m2 

In order to study the thermodynamical properties of the model it is interesting to 
observe that the Hamiltonian is left invariant by the action of a group of transformations. 
The group is described by G = Z2 x Z2 x Z2. 

We can represent a point in our parameter space as (m, J, h, a), where 

(mi \ I Jn J12 \ , I hi \ ^ I a 

m2 J \Ji2 J22 J \h2 I y 1 - a 

Therefore, given the limitations on the values of our parameters, the whole parameter 
space is given hy S = [-1, 1]^ x x R+ x x [0, 1]. 

If we consider the representation of G given by the 8 matrices 

ei \ ( ^ Vi\ 

, ej = +1 or - 1 and , 77, = +1 or - 1 

€2 y \ y 

we can consider the action of G on 5" as given by 

(i):GxS^S 

where 

0[M, (m,J,h,d)] = (Mm, MJM"^ Mh, MdM"^) 
for every x e S and M e G and it's straightforward to check that 

H(x) ^ H{cf>{M, x)). 

This can be easily done by writing the Hamiltonian per particle in vector notation as 

if(m,J,h,a) 1 /u - \ 

— ^ — — =— -(am. Jam) — (h,am). 

In order to obtain the analytic solution of the proposed model we consider the Boltzmann- 
Gibbs measure of weight 



and in particular we want to compute the average total magnetization per particle on 
that state 

For that purpose, it is useful to compute the pressure 



P = — logV 



One can show (see [13]) that in the thermodynamical limit (A^ oo) the pressure can 
be expressed as: 



P= sup/(/ii,/i2), (2.2) 



where 



/(/ii,/i2) = ^(^JiiaVi + -^22(1 - a)V2 + 2Ji2a(l - a)Ati/i2) + 
+ a/ii/ii + (1 - a;)/i2/i2 + 

,„(_i±-,„(i±-)_i^,„(i^)), 

+ (l-o)(-l±J!^\J l±^\ - 1.1 ( i^'l I . (2.3) 



The first two lines represent the internal energy contribution and the second the entropy 
in a state of magnetization = (/ii, 112)- 

Once we have the pressure it's easy to show that (m), in the thermodynamical limit, 
can be written as 

(m) = a(mi) + (1 - a)(m2) (2.4) 
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where (mi) and {1112) (the average magnetizations within the subsets Ji and I2) are found 
to be the maximizers of /(/ii,/i2) in (12.21) . 

The stationarity condition for the function /(/ii,/i2) in (12.21) gives 



= tanh(Jiia/ii + Ji2(l — a)/U2 + ^1) 

(2.5) 

/i2 = tanh(Ji2a/ii + J22(l - a)/U2 + /i2) 

This system has generically nine solutions, four of which are stable solutions corre- 
sponding to relative maxima. These can be found numerically by interpreting the (12.51) 
as a fixed point equation of a two dimensional map. Similar physical systems do appear 
in the litterature in the study of metamagnets (see [15j and [H]) but they only deal with 
a simplyfied subcase of the (12.51) which correspond to search the solutions on a subman- 
ifold of the entire phase space. Moreover the way the magnetic fields are coupled to the 
pre-interaction magnetizations is completely unusual and unprecedented in physics. See 
also ^16j for the use of the mean field ferromagnetic equations for a single population. 

In our model the values of hi and /12 shall not be considered as independent parameters 
but as functions of the average magnetizations and internal interactions in each original 
population when there is no mutual interaction between the two. 

Denoting by and the magnetization values at equilibrium within the population 
1 and 2 respectively, one has: 



3 Numerical results 

We have analyzed the numerical solutions of the system (12.51) and studied the behaviour 
in terms of the free parameters. 




(2.6) 



So our main quantity (m) is a function of five parameters: 



(m) = {m){a,Ju,J22,Ji2,ml,ml) . 



(2.7) 
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J-1 1 = 0.9 , J-22 = 1 .01 , m*-1 = -0.7 , m*-2 = 0.5 , alpha = fraction of immigrants 




Figure 1: (m) as a function of Ji^2 one? a 

The main quantity we study is (m). In the sociological context it may represent the 
average opinion of the interacting system consisting of the population of immigrants and 
residents. 

In particular, two main questions are addressed. First, what are the conditions that 
may lead to an instability of (m) in terms of the fraction of immigrants a? Second, how 
does the critical behaviour depends on the free parameters? 

As discussed in the previous section we may restrict our study to the case of Ji^2 > 0. 
The behaviour in the other regime can be deduced by symmetry. 

The results can be summarised in the following way. 

Figure [1] shows the total average magnetization (m) as a function of a and the mutual 
interaction strength Ji 2- 

What we see is that when J12 is small enough the magnetization is smoothly varying in 
a from 1712 to m* i.e. the value of the magnetization within the residents and immigrants 
before the mutual interaction takes place. But when the Ji,2 crosses a critical value we 
observe that the magnetization exhibits a discontinuous transition. The value of a at 
which the discontinuity occurs does not depend on Ji 2. This means that the critical Oc 
depends on the four parameters Ji^i, J2,2, ""^25 "^i- Numerical results show that ac can be 
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arbitrarily small for a suitable choice of the parameters it depends on. It is interesting 
then to investigate how ac depends on the interactions and on the original values of the 
magnetizations. 

We first present the results concerning the dependence of (m) in terms of the original 
magnetization within populations 1 and 2, i.e. the opinion of the immigrants and residents 
before interaction. In order to do so we study the following interaction matrix for j G [0, 7] 



(3.8) 



where Jn = J22 = 1 and J12 = j > 0. We work with a G [0, 1], and study (m) w.r.t. 
m* = {m\, ml^). 

The dependence of ac in terms of the cultural legacies m* appears to be rather intu- 
itive: the more the immigrant population is polarised (large negative values of m]^), the 
less is the amount of immigrants necessary to induce a phase transition. Equivalently, the 
more the resident population is polarised (large positive value of m\ ) the larger is a^- 

Pictures (a)-(f) in Figure [2] show some cases of surface Mm* (j, a) for m\ = 0.3 and 
several ml. The values of ac increase from to 1 when ml varies from -0.05 to -0.95 
(pictures (a) to (e)), while there is no abrupt transition when both ml and m"^ are positive 
(picture (f), top). The surface is symmetrical when taking both ml and negative 
(picture (f), bottom). 

Socially speaking, the transition has to be read as follows: for large enough intercul- 
tural interactions j, when a < ac the resident culture prevails, while when a > ac the 
immigrant culture prevails. An abrupt switch occurs when a = ac- Consequently, the 
critical value shows what fraction of immigrants is necessary to make the resident opinion 
lose its leadership over the entire population. 

The value of ac varies with m*. We may build a surface Sa^ (Figure [3]) where ac is a 
function of the non-interacting configuration m* with ml G [—1,0] and m2 G [0,1]. We 
may notice that etc = | when ml = — and the surface has the symmetry: ac{ml, 7712) = 
1 — ac{—m2, —ml). The value of ac is almost constant ac = when m* G Ti = { — 1 < 
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'Iast.daf using 1 :2:3 





Figure 2: Magnetization surface when (a) m* = (—0.05,0.3); (b) m* = (—0.2,0.3); 
(c) ui* = (-0.3,0.3); (d) ui* = (-0.5,0.3); (e) m* = (-0.95,0.3); (f) bottom graph: 
m* = (-0.3,-0.9); top graph: m* = (0.3,0.9). 



10 



'transiz.dat' using 2:1 :3 
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Figure 3: The transition value ac for the ferromagnetic model as a function of the non- 
interacting configuration (m^,m2) with m\ G [—1,0] and 777-2 ^ [O5 

m\ < — ^ , < 7Ti2 < — I — m\}. Due to the symmetry, = 1 when m* G T2 = { — ^ < 
m* < , I < 777^ < \-ml]. 

The a-critical surface Sa^ shows what is the amount of immigrants necessary to have 
a change, given the initial cultures of the two non-interacting populations. Its behaviour 
agrees with intuition: the stronger the immigrant original culture, the smaller the percent- 
age of them necessary to lead the opinion i.e. the critical percentage decreases with the 
immigrant culture strength. Viceversa the percentage increases with a stronger resident 
culture. 

We shall now investigate the dependence of ac on Ji 1 and J2,2- To do this we start 
from Figure [Hand increase Ji^i from 0.9 to 1.05, leaving the other parameters unchanged. 



The result of this variation is Figure 4(b) comparing it with the original picture 



shown on the side, we see how the discontinuity in the surface has been drawn forward to 
a bigger value of ac- the critical percentage turns out to be increasing with the strength 
of the immigrant interaction Ji^i. The more the immigrants tend to imitate each other, 
the less effective their influence is in the interacting system. This result may seem at odds 
with intuition, and suggests that a more diversified population is more likely to impose 
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= 0.9 , J-22 = 1 .01 , m'-l = -0.7 , m*-2 = 0.5 , alpha = fraction oT immigranls 

'1.dat' using 1:23 - 




J-1 1 = 1 .05 , J-22 = 1 .01 , m*-1 = -0.7 , m*-2 = 0.5 , alpha = fraction of immigranls 

'2.daf using 123 



--^ 0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 



alpfia 



(a) 




alpha 



m*-imm = -0.7 , m*-res = 0.5 , alpha - fractian of imm^rants 



(b) 




--^ 0.8 
- - 0.6 
^ 0.4 
0.2 


-0.2 
-0.4 
-0.6 
-0.8 



(c) 

Figure 4: Magnetization surfaces (a) before and (b) after increasing the immigrant popu- 
lation's interaction parameter Ju. (c) ac strictly increases with Ju. 
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m*-1 = -0.8 , m*-2 = 0.8 , alpha = fraction of immigrants 




Figure 5: Oc as a function of Ji^i and J2,2 



its cultural traits. 



we 



This turns out to be a general feature of our simulations. Consider Figure 4(c 
have set Ji,2 = 2 and we have let Jn vary between and 1.2. As a result, we see how etc 
increases monotonically with Jn. 

Similarly, by symmetry considerations, it is straightforward that an increase in the 
interaction strength within population 2 (J22) will lead to a decrease in a^- 

As a consequence, our study of the model reveals that an increase in the interaction 
strength within a population doesn't reinforce the population's own position within the 
total magnetization but, on the contrary, hastens the discontinuous transition towards 
the competing population's culture. The dependence of etc on Ju and J22 is summarised 
by Figure O 



4 Comments 

In this paper we have analyzed a two population mean field statistical mechanics spin 
model. The study of its phase structure in terms of the relative proportions of the two 
populations has been carried on. We have seen that the average magnetization of the two 
interacting populations may vary smoothly when the interchange coefficient is small, but 



13 



also abruptly when the coefficient is large. The critical value of the relative proportion 
has been studied in terms of both the initial values of the magnetizations within each 
population and a standard behaviour has been found. More intriguing is the observed 
dependence of the critical percentage in terms of the internal cohesion of each group. 
Due to a fine balance between internal energy and entropy we find that a strong cohesion 
penalizes the group. 

The model proposed is the simplest statistical mechanics model for the phenomenon 
of cultural contact, especially in the case of the residents and immigrants interaction. The 
dictionary associates a dichotomic opinion (like being in favour or against death penalty) 
to the two values of the spins. The interaction between two individuals is mapped into a 
ferromagnetic term in the Hamiltonian. The mean values within residents and immigrants 
of the country opinions are considered as cultural legacies and are compared to the average 
of the interacting mixed populations. The main result of this work from the social science 
perspective is to show that statistical mechanics predicts the possibility to have cultural 
dramatic changes during social contact. Moreover, it shows that the resident culture 
is more stable in its ability to survive the immigrant infiuence when imitation between 
residents is low and diversification is high. Equivalently, the power of the immigrant 
culture to take over and spread into the new country is lowered by a high internal imitation 
and low diversification. 

From the mathematics and physics point of view our model generalizes the theory of 
the metamagnets introduced in [13] and [T3]. The presence of extra degrees of freedom, 
in particular the pre-interaction magnetizations, in the two dimensional map makes our 
system difficult to study even though we may express the solution in a variational form 
like in the mean field theories: for instance it is not known yet what is the analytical 
condition needed to ensure the phase transition nor a detailed knowledge of its nature. 
Even simplified cases in which the solution is searched in submanifold of the entire phase 
space are not exactly solvable in a rigorous mathematical sense (see ^T^)] we will return 
on the rigorous results available in a future contribution [T3] . 
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